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Abstract 

The density matrix positivity is a natural counterpart of unitarity. The 
resulting constraints for various parton distribution and correlations are con¬ 
sidered. Their compatibility with leading order QCD evolution is guaranteed 
by the probabilistic interpretation of the latter. In the case of non-forward 
distributions the positivity constraints naturally imply the symmetric form of 
the evolution equation. 

The positivity of density matrix is just the positivity of its eigenvalues 

Ai > 0, (1) 

required, together with the normalization condition 

= ( 2 ) 

i 

by its standard probabilistic interpretation. As inequalities (Q) may be written in 
terms of quantities Trp* [Q 

Trp^ < 1 , 

2Trp^ — 3Trp^ > 1 , 

. ( 3 ) 

their preservation during interaction is guaranteed by the S— matrix unitarity, like 
the probability conservation (|^). Positivity is a natural counterpart of unitarity. It 
plays the more essential role, the more eigenvalues of density matrix are important. 

The positivity in QCD is pronounced in the framework of factorization, contain¬ 
ing the nonperturbative ingredient - parton distributions. They may be considered 
(at leading order) as a density matrices of partons in hadrons. The simplest appli¬ 
cation of positivity is resulting in the well known constraint for spin-dependent and 
spin-averaged distributions: 


f+{x),f_{x)>0, (4) 

or 

f+{x) - f-{x))\ < f{x){= U{x) + f_{x)) > 0. (5) 
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The non-diagonal (in helicity) elements of density matrix are also constrained, 
the well-known example being provided by the case of Soffer ineqnality for quark 
transversity distribution: 


\hi{x)\ < q+{x) = ^[g(x) + Ag(x)]. (6) 

The similar inequality 0 in the case of gluons is relating the contribution of 
twist 2 (G(a;)), 3 (AGr(x)) and 4:{Gl{x)): 

\AGt{x)\<^Ji/2G{x)Gl{x) . (7) 

It is most instructive to use this relation to estimate Gl (being the interesting new 
ingredient of nucleon structure from below: 

Gl{x) > 2[AGT{x)]yG{x). (8) 

This bound is analogous to well-known condition established long time ago by Doncel 
and de Rafael 0, written in the form 

IA2I A y/R , ( 9 ) 

where A 2 is the usual transverse asymmetry and R = a^/ax is the standard ratio 
in DIS. 

The similar bound relating the contributions of different twists may be derived 
for the twist-3 single spin asymmetries 0. Taking into account the so-called gluonic 
poles is resulting in the twist-3 term which is behaving like 1/(1 — x) with 
respect to the twist-2 unpolarized cross-section. This would necessarily mean the 
existence of twist-4 correction, behaving, in turn, like 1/(1 — As a result, the 
x—dependence of the asymmetry is governed by the generic equation [|. 

. 4 = ( 10 ) 

02 + 04/(1 — x)^ (1 — X)a2 + O4/ (1 — X) 

One can easily deduce two qualitative features of this result. 

i) The maximal asymmetry is achieved when the contributions to the unpolarized 
cross section of twists 2 and 4 are equal to each other, indicating the possible large 
contributions of higher twists. 

ii) This asymmetry (keeping twist 2 and 4 terms in the denominator) is twice 
smaller with respect to its ’’naive” estimate, when just the ratio of twist-3 and 
twist-2 terms is calculated. 

^In the case of single asymmetry of hadronic collisions, when the integration over momentum 
fraction of the unpolarized hadron is present, the average values of the coefficients a are assumed. 
If > 1 they are close to the values at x = 1 
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The same conclusions are valid, when Pt dependence of the asymmetry is con¬ 
sidered, 


b^M/pT &3 b^Mpx 

62 + b/^M'^ /p\ b2PT/M + b^M/pT PTb2 + M%4 


(11) 


where we recovered also the standard form of pt dependence. While maximal asym¬ 
metry is still estimated, up to a factor 2, by the ’’naive” asymmetry, its slope at 
small Pt governed by the twist-4 term. This is applicable for both fermionic and 
gluonic poles cases. Note that simultaneous change of x and px may be considered 
resulting in the same qualitative conclusions i) and ii). One may add to the posi¬ 
tivity properties of the observed asymmetries the bounds for the matrix elements, 
which are most simple for the gluonic pole case, when the matrix element T{x, x) is 
relevant: 


\T{x,x)\ < ^Jq{x)T{x,x, 0), (12) 

where T{x,x,0) is a twist 4 matrix element with 2 zero-momentum gluons. 

The positivity bounds may be derived i, I, 0, [Till also for the non-diagonal in 
momentum elements of parton-hadron density matrix, known as non-forward (off- 

0- 


forward) parton distributions |]12 


They are rather similar to 


although 


their exact form depends on the actual dehnition of such a distribution, existing in 
several versions. We present here the inequality for the quark off-forward distribu¬ 
tion 
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Iq{xi)q{x2) 

i-e 


and gluon non-forward distributions |ll3|. 


\xg{xuX2)\ < \JxiX2g{xi)g{x2), x=x-C- 


(13) 


(14) 


Here .^ = C(2 — C) is a ’’skewedness” parameter, while variables Xi ^2 measure the 
parton momentum fraction with respect to different hadron momenta. Their relation 
to the standard notation are different for off-forward {xi ^2 = / (1±.^)) and non¬ 

forward (xi = x,X 2 = {x — Q/{1 — ()) distributions. These variables, as was noted 


in 1^, are suitable for expressing the symmetry [14 resulting from T-invariance (and 
similar to the emerging in the case of forward twist-3 correlations P]), 


xg{x,x') = M{xi,X 2 ) = M{x 2 ,Xi), 


(15) 


(where M{xi,X 2 ) is the non-forward momentum distribution) and for t 
evolution: 


dM{xi,X2) _ P 

dt 271 Jxi 

M{xi,X 2 ) 
2 


—^^^—P{z,z)M{xi/z,X2/z') - 
^ dz T,, . P dz' /,,, 


0 1 — z 


/o 1 — z' 


InQ'^ 


(16) 


3 








































Here symmetric evolution kernel is related to the standard one p[|: P{z,z) = 
P{z',z) = z'{l — z)P{z,z), implying also the similar relation for diagonal ker¬ 
nel: P{z,z) = P{z) = z{l — z)P{z). Only the regular parts of the kernels are 
present here, while the singnlar contribntion is written down explicitly. To verify 
the symmetry of the evolntion eqnation one should note that new variables lead to 
the simple relation between the integrations over z and z' 


dz 


dz' 


Ixi Z{1 — z) Jx2 z'{l — z') 


(17) 


implied by the f-channel momentnm conservation 

Xi X2 1 1 

(1 — Xi)z (1 — X2)z' 1 — Xi 1 — X2 


(18) 


Snch a role of variables Xi,X 2 is not occasional, as they make non-forward distribn- 
tion look like X) -|- A*B (while forward ones are X H B'^) [0 and naturally 
manifest hermiticity, being the essential ingredient of the symmetry properties. The 
evolntion equation, in tnrn, should present this symmetry. It was stressed 0 
that in the case of the donble distribntion this symmetry is manifested, provided 
the overall factor 1 — (/2 [Q is extracted. However, the symmetry of nonforward 
distribntion is more complicated in that case. One can see that 


Ff(X) 


1-C 


c I 

1-C 


1-C 


)• 


(19) 


At the same time, the standard non-symmetric choice is leading to the cancella¬ 
tion of the overall factor 1/(1 — (C) in the r.h.s., so that resulting symmetry is just 
the interchange of Xi and X 2 -, as was shown above. 

The obtained evolution eqnation allows one to prove the stability of the positivity 
constraint (0) against Q^-evolution, following the general line of |^, 0. To do so, 
one shonld consider the positive qnantities (at some initial scale Qq) M±(xi,X 2 ) = 
aM{xi) + M{x 2 )/a ± 2 M{xi,X2), where a is an arbitrary positive nnmber. The 
ineqnality (|T^ 

\M{xi,X2)\ < \JM{xi)M{x2), (20) 


is just the result of minimization with respect to its variation, where M{x) = xg{x) is 
a diagonal momentnm distribntion. Making use of ([T^ , one may write the evolntion 
eqnations like: 


dM±{xi,X2) r dz 


dt 


2 tI Jxi z{l — z) 


{aP{z)M{xi/z) + P{z')M{x 2 lz')/a 


±2P{z,z)M{xi/z,X2/z')) - M±{xi,X2) 


"I dz 


1 -z 


mi- 


( 21 ) 


It is very important that the virtual contributions are diagonal in index ±, so that 
they cannot change the positivity of the distribntion (when the distribntion gets 
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too close to zero, it stops decreasing) like the exponentially decreasing positive 
function, which cannot change sign [^, [1^0. To prove positivity of the real term 
(c.f. p!6| ) it is sufficient to consider the minimization with respect to the variation 

of 2 ;-dependent (positive) a, which can only make the sum of two positive diagonal 
terms smaller, than in the actual case of minimization with respect to constant a: 


mm 


/ii + P{z')M{x 2 /z')/a ± 2P{z, z)M{x,/z, x^/z')) > 

!l^l^)^^^a{z){a{z)P{z)M{xi/z) + ^ 2P{z,z')M{xi/z,X2/z')) = 

^^{^P{z)M{x,/z)P{z')M{x2/z')±P{z,z)M{x^/z,x,/z')). ( 22 ) 


Writing down (^OD for xi —>• xi/z,X 2 —> X 2 Iz\ the sufficient condition of positivity 
of ( p2[) can be easily found 0: 


\P(z,z')\<^P(z)P(z'). 


(23) 


Such inequality is really valid 0 for all z, z', completing the proof of positivity. 
Although we considered here a pure gluodynamics, the mixing is improving the 


situation with positivity, providing extra positive terms, like in a forward case |T5 
One may consider the positivity of Q^-evolution in an ’’opposite” manner, when 
the general positivity constraints for various parton distributions may be used for 
obtaining bounds for yet unknown evolution kernels. From this point of view, the 
bounds (||,|D are most interesting, putting the constraints for twist-4 evolution. It 
is interesting, that the evolution of gri^x) = gi{x) -1- g 2 {x), proportional to A 2 , com¬ 
bines the twist-2 (from gi and Wandzura-Wilczek part of 5 ^ 2 ) and twist-3 terms, the 
latter, in turn, being multiplicative only in the limits a: —> 1 and Nc —» 00 jlQ 


However, the collinear singularities of real contribution to gT (essential for positiv¬ 
ity properties) after making use of symmetry due to T-invariance, are resulting in 
multiplicative renormalization. The virtual corrections should be cancelled in the 
non-diagonal in ± terms in evolution of grix) ± af(x) + fL{x)/a (c.f. 0 ))- 

The gluonic bound (^ may be used in a similar way to estimate the twist- 
4 evolution. Note that conclusion [Q about the absence of standard evolution for 
longitudinal gluon distribution is, generally speaking, applicable only for the effective 
twist-2 distribution, obtained by the contraction with the box gluon-photon diagram, 
the pole of which is cancelling the mass parameter, appearing, as usual, in the twist 
4 contribution. 

The account for NLO effects is more peculiar. The preservation of positivity 
for the parton distribution is depending on the choice of the factorization scheme 
and may be considered as an extra constraint for the latter [|T^. At the same time. 


full NLO result, dehning observable asymmetry, should respect positivity. Note, 


^Such property of virtual correction is making the positivity preservation especially simple when 
the evolution in x space is considered, while in the N space it requires the more elaborate analysis 
[p^, and the preservation of, say, Soffer inequality is conspired Q. 
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however, that the attributing the positivity in QCD to the physical asymmetry 
2m is, generally speaking, too strong condition. Say, Softer inequality is not re¬ 


lated to any physical asymmetry. Instead, it may be obtained from the positivity 
in the ’’gedanken-experiment” of the DIS scattering mediated by the hypothetical 
currents [j^. The one-loop coefficient functions of quark contributions (containing 
both g+(a:), hi(a;)) to this process, combined with the two-loop evolution kernels, 
should guarantee the validity of Soffer inequality at NLO. This would define the 
’’gedanken-experiment” factorization scheme, where the coefficient function of this 
process preserve its Born value and Soffer inequality is valid for the distributions. 
As to the ”Drell-Yan” factorization scheme, it should definitely preserve more weak 
inequality \h[x)\ < q{x). 

The inequalities 


may be used also for estimation of renormalon ^1[ contri¬ 
bution to twist-3 terms, the direct application being impossible because the Born 
term in the perturbative seriesis absent. Although positivity bound for A 2 in DIS 
is far from saturation, this is not necessarily the case for the twist 3 terms in single 
spin asymmetries, and the obtained bound may be of interest. 

To conclude, the positivity constraints may be used for both perturbative and 
nonperturbative ingredients of QCD factorization. They should provide the bounds 
for parton distributions and correlations, and for yet unknown evolution kernels. 


I am indebted to D.Boer, C. Bourrely, E. Leader, B. Pire, J. Soffer and R. Tanger- 
man for collaboration in the studies of described problems, numerous discussions, 
assistance and advice. 
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